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ABSTRACT 

' Context. The Galactic anomalous microwave emission detected between 10 and 90 GHz is a major foreground to CMB fluctuations. Well 
correlated with dust emission at 100 fjm, the anomalous foreground is interstellar but its origin is still debated. Possible carriers for this 
, emission are spinning, small dust grains carrying a permanent electric dipole. 

' Aims. To probe the origin of the anomalous foreground, we compare microwave data to dust IR emission on an angular scale of 1°, and search 
for specific signatures predicted by models of spinning dust. 

Methods. For the anomalous foreground, we use the 23 GHz all-sky map deduced from WMAP data by Miville-Deschenes and collaborators. 
' The infrared dust emission is traced by IRAS data. Models show that spinning dust emission is little sensitive to the intensity of the radiation 
, field (Go) for 10 < v < 30 GHz, while the mid-lR emission produced by the same small dust grains is proportional to Gq. To test this behaviour 
' in our comparison, we derive Go from the dust temperature maps of Schlegel and collaborators. 

Results. From all-sky maps, we show that the anomalous foreground is more strongly correlated with the emission of small grains (at 12 yum) 
, than with that of large grains (at 100 /im). In addition, we show that the former correlation is significantly improved when the 12 /jm flux is 
' divided by Go, as predicted by current models of spinning dust. The results apply to angular scales greater than 1°. Finally, from a model fit 
[ of the anomalous foreground, we deduce physical properties for PAHs that are in good agreement with those deduced from mid-IR spectroscopy. 
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1. Introduction 

As part of an effort towards accurate measurements of CMB 
fluctuations, experiments have motivated a detailed s t udy o f 



the Galactic foregrounds in the GHz-range. lKogtitetalJ ( ll996l) . 



iLeitch et all ( Il997h . and Ide Oliveira-Costa et alj ( 1 19971) found 



an unexpected emission excess between 10 and 90 GHz, 
which is correlated with dust far-IR but not with synchrotron 
emission. To avoid an inaccurate interpretation, this excess 
has been referred to as an anomalous foreground. Low fre- 



for V < 30 GHz Jde Oliveira-Costa et al. 


1999 


; Bandav et al. 


2003 ; Finkbeiner et alJ 2004 ; Davies et al. 


2006 


). Both this be- 



haviour and also the flux level of the excess are incompati- 
ble with what is known about the usual Galactic components 
in this sp ectral range: synchrotron, free - free, and thermal dust 
emission jde Oliyeira-Costa et al.lll999l 120041: lEagachd 12003 



Finkbeiner et alj 12004 iMiviUe-Deschenes et alj|2008l) . If the 
anomalous foreground is caused by spinning, small grains as 
proposed by Draine & Lazarian (1998, hereafter DL98) it 
should corre late more strongly with the m id-IR emission of 
small grains jde Oliveira-Costa et al ] l2002h than with the large 
grain (BG) far-IR emission. Until now, there is only incom- 
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plete evidence that this is the case, the difficulty being the sub- 
traction of zodiacal light in mid-IR data. Fr om a comparison 
of WMAP data to HI data, iLagachd j2003h showed that the 
anomalous foreground rises with decreasing column density, 
in a similar way to the emission of s mall, transiently heated 



j2006h 



grains. Towards the dark cloud L1622. lCasassus et al. 
found a stronger spatial correlation between the CBI 3 1 GHz 
and the IRAS 12 and 25 /dm bands than with the 60 and 100 
jum bands. In addition, models predict that the spinning dust 
emission is little sensitive to the intensity of the radiation field 
for 10 < V < 30 GHz (DL98, AU- Haimoud et al. 2 ( 309, Y sard & 
Verstraete 2009 hereafter YV09'). ICasassus et al.l j2008h tested 
this result in the p Oph molecular cloud: by comparing mid- 
IR (IRS, 13.3 yum) and 31 GHz (CBI) data, they showed that 
they are well correlated if the former is corrected for the excit- 
ing radiation field. In this paper, we further test the above pre- 
dictions of spinning dust models by comparing, on a 1 degree 
scale, enhanced IRAS data (IRIS) and the anomalous all-sky 
map of Miville-Deschenes et al. (2008, hereafter MD08). We 
also derive some physical properties of dust and compare them 
to those deduced from IR spectroscopy. 

The paper is organized as follows. Section 2 describes the ob- 
servational predictions of spinning dust models and how they 
can be used to probe the origin of the anomalous foreground. 
Section 3 presents the data sets used to reach this goal. Section 
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4 shows how the anomalous foreground correlates with dust 
emission. Section 5 lists some remarkable fields and presents 
the type of information that we expect to derive from the 
anomalous foreground study. Finally, Section 6 presents our 
conclusions. 



2. Behaviour of the spinning dust emission 

Nanometric-sized grains or Polycyclic Aromatic Hydrocarbons 
(PAHs) emit mostly in the mid-IR, whereas large grains dom- 
inate the FIR emission. The PAH emission is known to scale 
with the intensity of the radiation field GcQ dSellgren et al 



19851) . This is also true for the emission of large grains and for 
Go < 100. We show in Fig.[T]the behaviour of the spinning dust 
emission with Go in photometric bands predicted by our model 
(YV09). Models predict that the spinning dust emission near 
23 GHz is almost independent of Go when 0.01 < Go < 100 
(Ali-Haimoud et al. 2009, YV09, Fig. [B- This has strong ob- 
servational consequences. If anomalous foreground is spinning 
dust emission, we expect there to be a stronger correlation be- 
tween anomalous foreground at 23 GHz and IR emission di- 
vided by Go than with IR emission alone. Moreover this corre- 
lation should be stronger for 12 fim than 100 fim IRAS bands 
because the former traces the emission of small grains. 




1000 10000 



Fig. 1. Results obtained with the model of YV09, for a MRN 
size distribution (Nc - 30-216), where Nc is the number of C- 
atoms in the PAH and we assume a proton density of «// = 30 
cm \ Other gas phase parameters (temperature, and electron, 
proton, and ab undances) were computed with CLOUDY 
05.07 (IFerland et al. 1998) in the optically thin limit and as- 
suming 130 (320) ppm of carbon (oxygen) in the gas phase. 
Integrated intensity in the WMAP bands vs Gq: stars (red line) 
show the 23 GHz band, circles the 33 GHz band, squares the 
41 GHz band, diamonds the 61 GHz band, and triangles the 94 
GHz band. 



3. Data sets 

To carry out these correlations, we need maps of the anoma- 
lous and dust emission as well as for the Go values. IRAS 
is a natural data set to study dust IR emission. Our IR tem- 
plate is the new generation of IRA S images, called IRIS 
(IMiville-Deschenes & Lagachell2005b . in the 12 and 100 fim 
bands, which is corrected for most of the remaining instrumen- 
tal problems of the IRAS/ISSA dataset. Point sources were 
removed in the IRIS plates (at 5 arcmin resolution) usin g 
the method described in lMiville-Deschenes & Lagachel ( l2005h . 
The plates were then projected onto the Healpix grid, where an 
ecliptic-oriented filtering was applied to remove residual zodi- 
acal light emission (Miville-Deschenes et al. in preparation). 
Finally, the IRIS all-sky maps were convolved with a 1 degree 
FWHM Gaussian, smoothing out any imperfections related to 
the point source subtract i on. 

Miville-Deschenes et al. (l2008h performed a separation of 
components in the WMAP bands, using a physical approach 
to describe the Galactic foregrounds. We use this anomalous 
template at 23 GHz, inferred from their "Model 4". The main 
assumption made to obtain this map is that polarized emission 
at 23 GHz is dominated by synchrotron (no assumption about 
any correlation with dust). 

Fi nally, the Gp-m a p is de duced from the BG temperature map 



of iSchlegel et al.l (119981) that is inferred from the 140/240 
DIRBE ratio. We assume that the interstellar radiation field 
has the same spec tral distribution as the standard field of 
Mathis et al.l (Il983b . eve rywhere in the Galax y, and that the BG 
spectral index is /3 - 2 jDraine & Leelll984l) . The energy bal- 
ance of a single grain of size 0.1 fim then yields Go-{Tbc/^T -5 
K'f'^'^. All of these maps have been smoothed to the same an- 
gular resolution of 1°. 

4. Correlations 

Figure |2] shows the all-sky correlation of the 23 GHz anoma- 
lous flux with the dust IR emission. The anomalous foreground 
clearly correlates more strongly with the 12 //m band than the 
100 i^m (the Pearson correlation factor P is 0.90 and 0.82, re- 
spectively). A similar result was obtained by Casassus et al. 
(2006) towards the LDN 1622 cloud, but here it is the first time 
that it has been shown to also be true for the entire sky, follow- 
ing the removal of zodiacal light residuals at 12 jim. The corre- 
lation is also improved significantly when the dust IR emission 
is divided by Go (P - 0.90 to 0.95, in the case of the 12 fim 
bancd). This improvement concerns ~ 60% of the sky at 12 
fim. These regions are 1.4 to 1.6 times brighter at 23 GHz than 
the regions for which the division does not improve the cor- 
relation. However, in most of the regions where the division 
by Go does not improve the correlation, it also does not make 
it poorer It does only for 5% of the sky, which could be ex- 
plained by the uncertainties in the Go values. These correlations 
show the independence of the anomalous foreground of Go at 



' Scaling factor for the radiation field integrated between 6 and 13.6 
eV. The standard radiation fiel d corresponds to Gp = 1 and to an in- 
tensity of 1.6 X 10"^ erg/s/cm^ jParravano et alj2003h . 



^ Using a Monte Carlo method to simulate thermal noise in the 12 
fim map, we find that the Pearson coefficent 0.95 differs significantly 
from 0.90 with a confidence level greater than 99.9% (using a map 
containing 786432 pixels). 
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23 GHz and its link with the smallest grains. However, since 
the all-sky correlation is almost as strong with BG emission as 
with small grains, we are unable to draw firm conclusions at 
this stage. Across the entire sky, the emissions of PAHs and 
BGs are known to be correlated well. This is no longer true for 
particular fields, as we now discuss. 
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Fig. 2. All-sky correlations. Left panels show the anomalous 
foreground in the 23 GHz WMAP band versus the dust IR 
emission in the 12 and 100 yum IRAS bands. In the right panels, 
the dust IR emission has been divided by Go, the intensity of 
the radiation field. The contrast in point density between yel- 
low and black area is a factor 3. We indicate in each panel the 
value of the Pearson correlation factor P. 
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Fig. 4. Field centred on / = 343.6 + 5.25° and b = 22.4 ± 
2.6°. Correlations between a) anomalous foreground at 23 GHz 
and dust 12 fim emission, and b) anomalous foreground at 23 
GHz and dust 12 /im emission divided by Go. The black crosses 
show the low-Go area (median value equal to 1 .4) and the blue 
ones the high-Go area (median equal to 2.3). P is the Pearson's 
correlation factor 



5. Selected fields 



Anomalous emission, 23 GHz 




Fig. 3. Field centred on / = 343.6 + 5.25° and b = 22.4 ± 2.6°. 
The 23 GHz, 12 and 100 yum maps are in MJy/sr. 
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Fig. 5. A representative anomalous spectrum of our selected 
fields (taken from model 4 of MD08). The solid line is the 
fit with our model (YV09) comprising contributions from the 
CNM (10%, dashed line) and from the WNM (90%, dotted 
line). PAH parameters (see text) are m =0.3 D, and N^in = 24 
and 48 (CNM and WNM respectively). 

To test the spinning dust hypothesis further, we searched for 
fields of a few squared degrees according to the following cri- 



4 



N. Ysard et al.: Probing the origin of the microwave anomalous foreground 



teria: location outside of the Galactic plane, and bright at both 
23 GHz and 12 jim with Go variations as large as possible. 
Searching the sky maps by areas of 5° squared, we identified 
27 such fields. Figure |3] is an example of one of them. The 
anomalous and dust brightness maps correlate far more tightly 
when the latter is divided by gJI. The correlation plots indeed 
clearly illustrate two cases (Fig.H^) corresponding to different 
values of Go. The difference disappears when the 12 fim bright- 
ness is divided by Go (Fig. |4}5), as expected if the anomalous 
foreground is produced by the emission of spinning PAHs. 
For 5 of the 27 selected fields, we observe significant spatial 
variations between the 12 and 100 fim brightness maps (as 
shown in Fig. [3]). In these fields, we note that the correlation 
between the anomalous foreground and the BG 100 jum/Go is 
wors43 (Pearson's correlation factor P - 0.7 for the field in 
Fig. [3]) than with the smaller grains 12 //m/Go {P = 0.86). 
This shows that the anomalous foreground is correlated well 
with BG emission, only if BG emission is well correlated with 
IR emission characteristics of smaller grains. These results are 
consistent with spinning dust emission. 

We further test the spinning hypothesis and attempt to con- 
strain the electric dipole moment of PAHs in the selected fields. 
As discussed by YV09, the brightness of spinning PAHs at 23 
GHz, ^23, is given by Nh S pah m~€s, where Nh is the pro- 
ton column density, S pah is the abundance of PAHs, m is a 
scaling factor inferring the electric dipole moment of PAHs, 
fi(D) - m X N^^~ (where is the number of atoms in the 
PAH), and es is the rotational luminosity per solid angle and 
per PAH molecule. The PAH IR brightness in the 12 fim band, 
/i2, is proportional to Nh S pah Goej, where e/ is the IR lumi- 
nosity per solid angle and per PAH. The correlation coefficient 
between anomalous and IR brightness divided by Go is then 
m^xes/ei, where es depends on the number of carbon atoms in 
the smallest PAH molecules {N,„i„) and the fractions of neutral 
cold (CNM) and warm (WNM) diffuse gas. In Fig.|5] we show 
a representative fit to the observed anomalous foreground with 
our model. From the 27 selected fields, we find a mean ratio 
SiiKhilGo) - 1.3x10"^ withastandard deviation of4xl0"-'. 
Our model fits yield m -03- 0.4 D, N„,i„ =20-60 and about 
10% of CNM to account for both the 23 GHz and 12 jum emis- 
sion. These sizes are current ly invoked to explain the 3.3 Aim 



profil e in interstellar clouds (iVerstraete et alj|2001l : iPech et al 



2002h and the m-value is in good agreement with laboratory 
measurements for organic molecules (DL98). Thus, the rota- 
tional and vibrational emission of PAHs, as in current models, 
can consistently explain the anomalous and 12 //m emission for 
plausible properties of PAHs. 

6. Conclusions 

From an all-sky, degree-scale comparison of the 23 GHz 
anomalous map with dust IR emission, we have found that the 



^ The improvement in the correlations is significant to a confidence 
level greater than 99.7% for the 27 regions. 

* The correlation is also poor with the 60 yum / Go and for 3 of 
them with 25 ^m / Go (P = 0.18 and P = 0.7 for the field in Fig.[3l 
respectively). 



anomalous foreground is well correlated with the 100 jum IRAS 
band. Using an enhanced set of IRAS maps, we have shown for 
the first time that the anomalous foreground is correlated with 
the 12 yum band across the entire sky and that the correlation 
is tighter than with the 100 flux. This correlation becomes 
even tighter when the 12 jum flux is corrected for the intensity 
of the radiation field Go, indicating that the anomalous emis- 
sion is independent of Go at 23 GHz on a I degree scale. These 
findings strongly argue in favor of a spinning dust origin to the 
anomalous foreground. Current models predict that the spin- 
ning dust emission is dominated by the smallest dust grains 
(PAHs) carrying the 12 yum flux and that the corresponding 23 
GHz emission is almost independent of Go. From a model fit 
of both microwave and IR data in selected fields with strong 
Go contrast, we deduce the physical properties of PAHs (sizes, 
electric dipole moment) that are in good agreement with results 
obtained from mid-IR spectroscopy. 
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